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Abstract 



We evaluate in perturbative QCD the semi-inclusive Drell-Yan cross-section for the production of a single 
hadron accompaining the lepton pair. We demonstrate to one loop level a collinear factorization formula 
within the fracture functions approach. We propose such a process as a factorization analyzer in hadronic 
collisions. Phenomenological implications at the hadron colliders are briefly discussed. 
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Introduction 



The Drell-Yan pQ process has unique features among high energy hadronic reactions. The mea- 
surament of the relative cross-sections represents in fact a basic normalization in the LHC physics 
program, in particular by proving the universality of parton distribution functions as measured in 
Deep Inelastic Scattering and by possibly pinning down non-standard evolution in the initial state 
parton cascade at small values of Bjorken variable x. 

Universality of parton distribution functions follows directly from the factorization theorem 
for the Drell-Yan process. Despite some initial controversy on such a factorization in the early 
'80, a series of papers [2] supported this hypothesis with specific model calculations and then a 
factorization proof was finally given in Refs. [31 H]. It was shown in particular that, when diagrams 
with soft gluon exchanges between active and spectators partons are summed over all possible cuts, 
soft gluon contributions decouple from the short distance cross-section, convalidating factorization 
at parton level as firstly conjectured by Drell and Yan. All the complications in proving the 
factorization can be related to the presence of two hadrons within the initial state. 

The Drell-Yan process has furthermore two more appreciable features with respect to other 
short-distance hadron-induced cross-sections: a) the perturbative scale can be accurately recon- 
structed by measuring the invariant mass Q 2 of the lepton pair and b) the final state is free from 
Q CD-corrections thus providing a clean tool for the study of initial state radiation pattern. 

Given these properties, it may become natural to investigate what happens if an additional 
hadron, accompaining the lepton pair, is identified in the final state: 

Pi + P 2 ^7* + h + X , (1) 

where Pi and P2 denote the incoming hadrons, h is the identified hadron in the final state, 7* is the 
virtual (eventually electroweak) boson and X the unobserved part of the final state. The idea of 
using the Drell-Yan process as a perturbative trigger were first investigated in Ref. [Sj and will be 
fully exploited also at the LHC [SI E] • If the invariant mass Q 2 of the lepton pair is large enough so 
that perturbation theory applies, following the arguments of Ref. [8] , the factorization property of 
the cross-section for the process in eq. ([I]), should depend on the region of phase space in which the 
final hadron h is detected. In particular if h is produced at sufficiently high transverse momentum, 
p 2 h±) then in such phase space regions (we refer to them as to central) standard perturbative QCD 
technique should be applicable, see for istance Ref. [SJ. If h has instead a low p\ ± and thus is 
detected in the so-called target fragmentation region, arguments against factorization have been 
already given in Refs. pH [HI [12] . 

It is therefore highly desiderable to have a standard perturbative framework which can be used 
as a "factorization analyzer" in both region of phase space. The first step is thus to provide a 
parton model formula which accounts for the production of an additional hadron. In the inclusive 
case the cross-section for the production of a lepton pair of mass Q 2 in the collision of two hadrons, 



1 



of momenta Pi and P2, can be written as [I] 



da DY (r) _ Ana 2 
dQ 2 ~ 9SQ 2 



xix 2 



T 



), (2) 



with S the hadronic center of mass energy, S = (Pi + P2) 2 , and r = Q 2 /S as in Ref. [32]. The 
sum runs over all quarks and antiquarks flavour but not on gluons which can not directly couple to 
electroweak bosons. The parton distribution functions f q (x) depend on the fractional momentum 
of the parton entering the hard scattering. In particular no scale dependence is indicated as 
appropriate in the naive parton model formula. In such an approach to the semi-inclusive Drell- 
Yan process of eq. flTJ, where QCD higher order corrections are absent, we assume that the final 
state hadron is "non-perturbatively" produced in the target fragmentation region of Pi (T^rJ or 
P2 (7^t 2 ) by means of a fracture function MJ l ,p(x, z). These distributions give the conditional 
probability of finding a parton i with a fractional momentum x while an hadron h, with fractional 
momentum z of the incoming hadron momentum P, is detected in the target fragmentation region 
of P, see Ref. [H]. The collinear and soft factorization of these distributions in semi- inclusive DIS 
has been proven respectively in Refs. [T5l[T2] and by an explicit 0(a s ) QCD calculation in Ref. [16J. 
Supported by these results, continous efforts have been devoted to the extraction diffractive parton 
distributions from HERA data. A combined analysis of both the diffractive and leading proton DIS 
data in terms of fracture functions has been presented in Ref. [17] while for more recent analysis 
we refer to Refs. [TS]. Since no perturbative emissions are allowed in a parton model formula, 
we assume that "bare" fracture functions describe hadron production in the target fragmentation 
region Ttr 1 of Pi if 9 cm = and in 1Zt 2 of P2, if cm = n, where 9 cm is the relative angle between h 
and Pi in the hadronic center of mass frame. Phase space separation between target and central 
region is indeed unphysical and in this particular case also frame dependent. However this choice 
can be shown to be the more suitable in order to prove the factorization of the collinear singularities 
of the cross-sections. In the following we will consider the next-to-simple differential cross-sections 
for producing a lepton pair of invariant mass Q 2 3> Aq CD , accompained by an additional hadron h 

with fractional energy z = 2Eh/y~S (defined in the hadronic center of mass frame) and integrated 
over its transverse momentum, p 2 h± . By defining the combination 



a straightforward generalization of eq. ([2]) leads to the parton model formula for the semi-inclusive 
Drell-Yan process: 



According to momentum conservation, the convolution integrals in eq. (j4]) must satisfy the con- 
straints 1 — z > X\%2 > t and 
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da DY {r) _ Ana 2 
dQ 2 dz ~ 9SQ 2 



X\ + z < 1 in !Z Tl and x 2 + z < 1 in 1Zt 2 > 
X2 + z < 1 in IZt-l and x\ + z < 1 in 1Zt 2 , 
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Figure 1: A pictorial rapresentation of eq. (JH). 



for the first and the second terms in square brackets of eq. (jlj) respectively. The number of terms 
appearing in eq. (j3J) is twice the number appearing in the inclusive case since each fracture func- 
tions selects its own fragmentation region. Eq. (j4]) represents for the moment only a factorization 
conjecture for the process we are considering and is sketched in Fig. (CQ). Our purpose is to show 
that eq. Q does survive to the inclusion of radiative corrections and that the factorization of 
collinear singularities is possible when the proper subtraction terms for bare fracture functions are 
added. The present QCD-based calculation deals with the standard soft gluon exchange between 
active partons but it is blind to soft gluon exchange between spectators whose effects are taken into 
account in Refs. [2j [3J S] and could spoil the factorization, as suggested in Ref . p2] . In the presence 
of such factorization breaking effects, as emerged in recent phenomenological analysis, diffractive 
parton distribution extracted from HERA data do overestimate diffractive cross-sections measured 
at Tevatron by an order of magnitude [19J. Fracture functions appearing in eq. (jlj) therefore can not 
be related to the ones extracted from Semi-Inclusive DIS data. In a strong factorization breaking 
scenario the factorized form M ® / itself could in principle be questioned, since the exchanges of 
low momemuntum gluons cannot be uniquely absorbed neither in the definition of fracture nor 
of parton distribution functions so that they pertain to the reaction as a whole. The formalism 
we use can not indeed penetrate inside the details of soft factorization, however it may constitute 
a quantitative next-to-leading order guideline for estimating the magnitude of its breaking when 
passing from expected non-factorizing phase space region at low p\± to expected factorizing ones 
at high p1 ± . The possible identification of an intermediate scale or range of scales at which this 
transition may occur would constitute an important matter of inspection insight in the dynamics 
of the factorization mechanism. 

Evaluation of NLO corrections 

As outlined above, we have adopted the normalization and conventions as of Ref. [13] . The cross- 
sections for the leading order partonic sub-process q(pi) + q{p2) - > 7* is defined by: 




(5) 
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with Nc the number of colours and w = Q 2 /s with s = (p\ + p^) 2 the partonic center of mass 
energy. The space-time dimension is n = 4 — 2e. Eq. (jSJ) is implicitely contained in parton model 
formula, eq. ([2]) and eq. (T4|). Moving to next order in perturbation theory requires the evaluation of 
the real emission diagrams q + q —>■ 7* + g and q + g —>■ 7* + g along with virtual corrections to the 
Born amplitude. The results for the matrix elements squared, averaged over initial state colours 
and spins and summed over final ones, as well as the virtual contribution, can be found in Ref. [13J. 
Both the real and virtual contributions depend on the renormalization scale p 2 . Moreover, the 
cross-section dependence on the variables associated to the produced hadron is entirely contained 
in fracture functions for both Born and virtual contributions. 

So far no new ingredients have been added to a standard NLO perturbative calculations, the 
major subtleties coming from the two-particle final state of real emissions diagrams. In order not 
to obscure the renormalization procedure we restrict ourselves to the discussion of the qq channel 
which indeed contains all the essential part of the calculations. In the following we will demonstrate 
that, at order 0(a s ), we can organize radiative corrections in such way that uncancelled (collinear) 
divergencies can be reabsorbed into the "bare" distributions appearing in eq. (jlj), in analogy with 
the semi- inclusive DIS calculation as performed in Ref. [16] . 

When the final state hadron is observed in IZ^ or 1Zt 2 (for istance in hard single diffract ive 
events in hadronic collisions), we assume that it has been produced non-perturbatively from fracture 
functions. All the perturbative real radiation thus must be integrated over and virtual corrections 
added. The corrections to the semi-inclusive Drell-Yan process in the target fragmentation regions 
turn out to be the same as in the inclusive case. The two body phase space can be easily integrated 
in the partonic center of mass frame. It can be shown by explicit calculations that the double poles 
in the e-expanded results cancel in the sum of real and virtual term. This latter feature is peculiar 
of the qq channel. For the singular contributions to the cross-section in the target regions we obtain 



daP Y (r) 4na 2 f 1 2 dx 1 f 1 dx 2 



dQ 2 dz 9SQ 2 J T X\ J j_ x 2 



e 2 q [m^xx, z)f q {x 2 ) + (x! <-> x 2 ) 



e 2tt qqy J \ Q 2 J r(l-2e) 



T(l-e) 



(6) 



with w = t/xix 2 and P%j{w) the Altarelli-Parisi splitting function. Next we consider the production 
of the observed hadron h by the fragmentation of a, real, final state parton (i.e. a gluon) in the 
partonic sub-process. This 0(a s ) production mechanism is sketched in Fig. (J2J). We address this 
phase space region as "central", where the hadron h is allowed to be produced at high p 2 ± . In 
the collinear limit the "central" region collapses to the target fragmentation one. In general such 
correction is expected of the form [20J 



da® Y (t) 4na 2 2 [ dp f dx\ f dx 2 t ( w , XT ^hi / \ dcr%, 
dQ 2 dz = 9SQ 2 



The variable p = 2Ej i /\/S is the partonic analogue of z, being Ek the energy of the outgoing 
gluon and da 9 qq - the differential partonic cross-sections in the considered channel. The phase space 
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Figure 2: Hadron production in the central region 



integrations are performed in the hadronic rather then in the partonic center of mass frame. The 
final state gluon depends on its fractional energy p and on an angular variable y — (1 — cos# cm )/2, 
being 8 cm the relative angle between k abd P\. These variables are however not independent and 
constrained by the formula 

/ x = 2(xix 2 - r) 

P{ ' x l + x 2 + (2y-l)(x 1 -x 2 ) ' 1 ' 

In order to keep compact the expressions it is useful to rewrite the p convolution as a y integral, 
with p depending on y via eq. ([5]). The central contribution to the cross-sections therefore reads 

<icr? y (r) Ana 2 y-^ 2 f 1 ^ f 1 dx\ f 1 dx 2 ^ ^ \]j h ( I * (9) 

dQ 2 dz 9SQ 2 q J J ri (r,z;y) %1 J r 2 {r,z;y) X 2 " 1 " 2 9 dx X dx 2 dy p' 

The integration limits ri and r 2 are obtained by imposing momentum conservation and the condition 
p > z in order to guarantee that the energy of the parent parton is greater than the energy of the 
observed hadron: 

/ n T + z(l-y) r + Xl zy 

n{T,z;y) = — , r 2 (r,z;y) = -. (10) 

l-zy xx-z(l-y) 

The emitted gluon in this region it is thus not allowed to be soft since it is required to produce the 
observed hadron. Once the invariants t and u appearing in the relevant matrix element squared 
are specified in this frame, we perform a standard e-expansion of the result in two disjoint singular 
limits, i.e. for y — > 0, 1. Retaining only 0{e~ 1 ) in the expansions we get the singular contributions 
to the cross-sections in this region of phase space: 

d(Tc Y {r) = Ana 2 C 1 dxi C 1 ^2,/ )f( )n h( zx 2 
dQ 2 dz 9SQ 2 J T+Z Xl J T/{x ^ z) x 2 Jq{Xl> MX2j Ax^-r) 2n 

( 1\ Hi -e) ( knjZ Yf, ( r \ x 2 
A e /r(l-2e) V ; r ^^\ Xl x 2 ) x lX2 -r 
Ana 2 f 1 dxi f 1 dx 2 ... . h ( zx\ \a s (p 2 ) 

f q {Xx) fq{x 2 )D n \ 



9SQ 2 Jr/(i-z) %i J z+T / Xl x 2 9 \x 1 x 2 -rJ 2n 

1\ r(l - e) ( Anp 2 r y p / r \ x x 
' e)V{\-2e)\ Q 2 ) m *~ q \ Xl x 2 ) x x x 2 - r ' 1 J 
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The first term in eq. (1111) is singular for y — > while the second for y 
do appear unregularized Altarelli-Parisi splitting functions, P( 



-> 1 . In the previous expression 
'w), since no interference with 



virtual contribution is present. Eqs. ([6]) and eq. (fTTT) thus represents all the singular collinear 
contributions to the semi-inclusive Drell-Yan cross-sections. If we were considering the inclusive 
Drell-Yan case, the only singularities would be the ones shown in eq. ([6]). The subtraction of 
singular term in the partonic cross-sections would be performed by lumping the divergence into the 
bare parton distributions /. In the MS subtraction scheme it reads: 



/<(0 



1 du 



u 



SijS(l — u) + 



W/i 2 ) r(l-e) 
e 2vr r(l-2e) 



/A 



£ 2 



(12) 



In the previous equation /i 2 and /i 2 are respectively the renormalization and factorization scale. 
Since any observable built using eq. (TT2"j) can not depend on fi 2 , the derivative with respect to 
ln/i 2 gives standard QCD evolution equations [21] • In the semi- inclusive case however there are 
additional singularities, eq. ffTTl) . and the subtraction of eq. ([12]) would not be sufficient to render 
the cross-sections infrared finite. Fracture functions, however, have been shown to have a more 
complex evolution equations with respect to the one of /. They contain an inhomogeneous term 
which accounts for hadron production in the target fragmentation region of the projectile by the 
showering of initial state radiation [13]. The analogous of eq. (|12|) for bare fracture functions is 
obtained in the context of semi-inclusive Deep Inelastic Scattering 



+ 



du 

e/(i-o u 

du 



l a s (tf) T(l-e) / 47r^ y 
U)+ e 2tt 1(1-26)1 y? ) nAU) 



MC 1 

L wla s (/i 2 ) r(l 



J \u 



+ 



(^K 

u l-u£e 2tt T(l-2e) V /x 2 



(fc)i^-j 



(13) 



The first term on r.h.s of eq. f)13p has the same subtraction structure as for parton distribution, 
eq. (1121) . The singularity is due do collinear radiation accompaining the active parton, while the 
hadron in the final state is non perturbatively produced by the fracture functions itself. In the 
second term of eq. (fl3"j) instead, the singularity is due to the observed hadron being collinear to the 
incoming hadron and generated by the perturbative fragmentation of the emitted collinear parton. 

At this point we insert eq. f|T2|) and eq. f|T3l) in the parton model result, eq. (j3]). After some 
algebra and integral manipulation is easy but lenghtly to show that eq. f|T2l) and the homogeneous 
term of eq. f|T3l) produce exactly the pole term in eq. (jSJ) with the opposite sign. The inhomogeneous 
term in eq. ([TBI instead reproduce the singular term in eq. ffTTl) . All the singularities therefore cancel 
and we are left with 
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da DY {r) 
dQ 2 dz 




The function C qq and K qq are infrared finite and depend explicitely on the factorization scale n F . 
All bare distributions are replaced by renormalized as indicated by the explicit factorization scale 
dependence. In particular C qq is the same as in the inclusive Drell-Yan case whereas K qq - is specific of 
the semi-inclusive process. Its explicit form along with the coefficient for the gluon initiated channel 
will be reported in a separate paper, as well as a more detailed description of the calculation. In 
eq. (|14p we may set fi F = Q 2 in order to remove potentially large logarithmic corrections of the 
type ln(fi F /Q 2 ) from the coefficient functions and resum them by using the appropriate evolution 
equations for / and M. 

Conclusions 

We have shown, with a fixed 0(a s ) perturbative QCD calculation, that the partonic cross-sections 
for the semi-inclusive Drell-Yan process do factorize at the collinear level. We therefore implicitely 
confirm the general widespread idea indicating that soft exchanges between active and spectators 
partons [10] as responsible for factorization breaking in semi-inclusive reactions. 

We would like to conclude by listing a few remarks and proposals. As opposed to full inclusive 
observables, semi-inclusive ones in hadronic collisions are affected by soft gluon exchange contribu- 
tions and could therefore act as factorization analyzers at phenomenological level. In the central 
region, high p\ ± hadron production should follow the pattern predicted by perturbative QCD. 
When the detected hadron is, instead, a low p\ ± proton, diffractive processes may occur [22J and 
the partonic structure of the color-singlet exchanged object may be studied [23] • In these partic- 
ular cases a non universality of diffractive parton distributions, as taken from diffractive DIS and 
hadronic collisions, was suggested in Ref. This was experimentally reported in Ref. [19J. It 
would be interesting to establish with a comparison with the data whether a factorization breaking 
shows up only in a diffractive kinematic regime or if it manifests itself also in processes with a gap- 
less final state containing a single hadron as well in the target fragmentation region. For the same 
reason it would be interesting also to study, within the proposed approach, light mesons production 
which is sensitive to the soft, high multiplicity, fragmentation process. 
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The present work can be also generalized to double hadron production. The evaluation of a dou- 
ble hadron production cross-section needs a full 0(a 2 s ) QCD calculation. However, an approximate 
result could be obtained if one considers two hadron at low p 2 h± observed in opposite fragmentation 
regions with respect to the incoming hadrons. In this case higher order corrections for this process 
should be the same as for inclusive Drell-Yan process, when the proper kinematics is taken into 
account. In Ref. [24] we have suggested an analogous formula, in leading logarithmic approxima- 
tion, for double-inclusive Drell-Yan production which also includes the additional dependence on 
the invariant momentum transfer at the proton's vertex t\ and £2- 

Finally we are thinking to a generalization of the present approach to include gluon initiated 
hard processes [25J whose relevance in diffractive Higgs production was first suggested in Ref. |26j. 
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